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Abstract 


Our  recent  paper  [6]  presented  a  model  for  thin  plates  with  rapidly 
varying  thickness,  distinguishing  between  thickness  variation  on  a  length 
scale  longer  than  ("a  <  1").  on  the  order  of  ("a  “  1"),  or  shorter  than 
("a  <  1")  the  mean  thickness.  We  review  the  model  here,  and  identify  the 
"a  <  1”  case  as  an  asymptotic  limit  of  the  case  "a  =  1".  We  then  present 
a  convergence  theorem  for  the  "a  -  1"  case,  showing  that  the  model  correctly 
represents  the  solution  of  the  equations  of  linear  elasticity  on  the 
three-dimensional  plate  domain,  asymptotically  as  the  mean  thickness 
tends  to  zero. 
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1.  Introduction 


In  (6]  we  presented  a  model  for  the  bending  of  symmetric,  linearly 
elastic  plates  with  rapidly  varying  thickness.  The  model  distinguishes 
between  three  cases,  in  which  the  thickness  varies  on  a  length  scale  longer 
than  ("a  <  1")*  on  the  order  of  ("a  -  1")*  or  shorter  than  ("a  >  1")  the 
mean  thickness.  (See  §2B  for  the  precise  meaning  of  the  parameter  a.) 

Our  main  goal  in  this  paper  is  a  convergence  theorem  for  the  "a  =  1" 
case,  assuming  that  the  thickness  varies  periodically  and  the  plate  edges 
are  clamped.  The  corresponding  convergence  result  for  flat,  homogeneous 
plates  is  well-known  [8,9].  Caillerie  has  recently  proved  an  analogous 
result  for  flat  plates  with  rapidly  varying  composition  [2,3];  his  method 
does  not,  however,  extend  readily  to  the  problem  treated  here. 

We  begin,  in  sections  2  and  3,  with  a  review  of  the  model.  The  vertical 
midplane  displacement  satisfies,  in  each  case,  a  fourth-order  equation 


3  q(M  ,  3  w)  =  F  ; 
a0  aByd  yd 


the  formula  for  the  effective  rigidity  tensor  M  „  .  differs,  however,  in 
the  three  cases.  Proposition  3.1  identifies  the  "a  <  1"  effective  rigidity 
as  an  asymptotic  limit  of  the  "a  **  1"  case.  It  is  possible,  at  least 
formally,  to  obtain  the  "a  >  1"  effective  rigidity  from  the  "a  ■  1"  case 
as  the  opposite  limit,  in  which  the  horizontal  scale  of  the  cell  be¬ 
comes  small.  To  make  this  analysis  rigorous  one  seems  to  need  regularity 
assumptions,  the  validity  of  which  are  not  obvious  in  any  interesting 
cases.  Though  this  analysis  is  not  included  in  the  present  paper,  we 
mention  it  here,  in  connection  with  Proposition  3.1,  as  an  indication  that 
the  case  ”a  *  1"  -  when  the  period  is  comparable  to  the  mean  thickness  - 


is  in  a  certain  weak  sense  universal. 


2 


T 


Sections  4-6  present  the  convergence  proof.  Our  principal  tools  are 
two  integral  inequalities,  an  "averaging  lemma",  and  lots  of  integration 
by  parts. 

The  first  integral  estimate.  Proposition  4.1,  is  a  Korn-type  inequality, 
in  which  the  dependence  of  the  "constant"  on  thickness  is  made  explicit. 

The  second  estimate.  Proposition  4.2,  asserts  a  weak  form  of  Kirchhoff's 
hypothesis  for  any  displacement  field  which  obeys  certain  symmetries. 

Our  "averaging  lemma".  Proposition  5.1,  amounts  to  an  estimate  for  the 
work  done  by  a  vertical,  locally  equilibrated  load.  We  call  it  an  averaging 
lemma  because  it  is  used  in  the  convergence  argument  to  replace  oscillatory 
expressions  by  their  local  averages,  modulo  error  terms. 

The  main  convergence  argument,  given  in  section  6,  proceeds  very 

straightforwardly:  we  simply  estimate  the  energy  of  u.  “  u  ,  where  u 

solves  the  equations  of  three-dimensional  linear  elasticity,  and  u*  is  an 

ansatz  motivated  by  formal  asymptotic  expansion.  The  main  result.  Theorem  6.1, 

is  thus  one  of  convergence  in  energy  on  an  e-dependent  domain.  The  conver- 

2 

gence  of  the  mean  vertical  displacement,  in  a  weighted  L  -norm  on  the  mid¬ 
plane,  follows  as  a  corollary. 

The  convergence  analysis  presented  here  applies  only  to  the  case  "a  =  1", 
and  only  to  plates  with  periodic  thickness  variation  and  clamped  edges.  The 
method  appears,  however,  to  be  rather  more  general.  We  believe  it  could  be 
applied  with  other  boundary  conditions  at  the  plate  edge,  and  with  plates 
whose  thickness  is  "locally  periodic"  or  "quasiperiodic"  in  the  sense  of 
{6].  An  analysis  of  the  cases  "a  <  1"  and  "a  >  1"  could  perhaps  be  done 
following  a  similar  outline. 


Structural  engineers  are  interested  in  plates  of  the  type  studied 
here,  because  they  may  be  stronger  per  unit  weight  than  uniform  or  slowly 


3 


varying  ones,  in  certain  design  contexts.  Some  references  to  the  literature 
on  structural  optimization  can  be  found  in  [6].  It  is  natural  to  ask  which 
scaling  -a<l,a*l,ora>l-  produces  the  strongest  structure:  we 
hope  to  address  this  issue  in  a  forthcoming  paper. 

We  acknowledge  with  pleasure  advice  from  George  Papanicolaou  on  aspects 


of  this  project. 
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2.  Preliminaries 

3 

We  shall  write  x  -  (xj.x^x^  for  vectors  in  R  and  x  -  (x^x^ 

2 

for  vectors  in  R  .  Integer  latin  indices  will  range  from  1  to  3,  and  Greek  ones 

from  1  to  2;  the  summation  convention  applies  whenever  indices  are  repeated. 

2 

We  write  3t  =  3/3xt  and  3^  -  3  /Sx^x^  . 

2A.  Constitutive  laws 

3 

Associated  with  any  displacement  u  -  (u^.^.u^)  of  R  is  its  linear 
strain  tensor 

(2.1)  =  ^8jui+8iuj> 

and  the  corresponding  stress  tensor 

(2. 2)  =  Bljk£eu(u)  . 

The  fourth-order  tensor  s8**8^1-®8 

®ijk£  ®j ik£  “  Bij£k  =  ®k£ij  ; 


we  assume  that  the  elastic  energy 

®ijk£  eij  *k l 


is  positive  definite  on  symmetric  tensors. 

We  shall  always  assume  that  the  horizontal  planes  are  planes  of 
elastic  symmetry;  this  means  [7} 


Ba8y3 


0 


Ba333 


0  . 


Finally,  we  define  a  positive  definite  fourth-order  tensor 
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(2.3) 


B  =  B  -  Ba333  BY»3 
cxByd  '  af5Y6  B3333 


2B.  Plate  geometry 

The  plate  geometry  is  determined  by 

(2.4a)  a  smoothly  bounded  domain  71  in  the  x^  -  Xj  plane, 

representing  the  midplane; 

(2.4b)  a  real  parameter  a,  0  <  a  <  00  ,  determining  the 

length  scale  of  the  thickness  variation,  and 

(2.4c)  a  bounded  function  h(n)  *  0  ,  defined  for  any 

2 

n  €  R  and  periodic  in  na  with  period  , 

a  =  1,2  . 

The  three-dimensional  region  occupied  by  the  plate  is 

R(e)  *  (x  :  x  i  JJ  ,  | x3 1  <  eh(x/ea)}  ; 

R(e)  denotes  its  natural  periodic  extension 

(2.5)  R(e)  ■  (x  :  x  €  R  ,  |x3|  <  eh(x/ea)}  . 

2  a 

We  assume  throughout  that  R(e)  is  a  connected,  C  ’  domain,  for  some 
Holder  exponent  a  >  0  .  The  function  h  may  nonetheless  have  discontinuities  - 
i. e.  parts  of  3R(e)  may  be  vertical;  and  h  may  vanish  on  a  set  of 
positive  measure  -  i.e.  our  plates  may  have  holes.  (In  section  3D,  where 
we  study  an  asymptotic  limit  of  the  a  ■  1  case,  we  shall  impose  additional 
smoothness  assumptions  on  h  .) 

We  denote  by  3qR(e)  the  outer  edge  of  the  plate. 
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8qR(e)  =  {x  :  x  €  30  ,  1x^1  <  eh(x/ea)  ; 

3+R(e)  and  3_R(e)  are  the  remaining  parts  of  3R(e)  above  and  below 
Q  ,  respectively;  vE  is  the  outward  unit  normal  to  3R(e)  . 

When,  in  the  following,  we  call  a  function  "periodic  in  n  "  we  shall 
always  mean  that  it  has  the  same  periods  L  =  (L^.L^)  as  h  .  It  will 
often  be  necessary  to  average  a  periodic  function  f(n)  with  respect  to 
n  : 


Mf 


L1L2 


We  shall  use  the  norm 


full. 


£ 


L1  fL2 

f(n)dn.dn-  . 

0  J0 


2 

The  same  notation  will  be  used  for  tensors,  in  which  case  |g|  denotes 
the  sum  of  the  squares  of  the  components. 

2C.  Loads  and  equations  of  equilibrium 

The  following  disucssion  applies  for  a  *>  1;  when  a  f  1,  it  is  more 

natural  to  work  with  the  load  per  unit  projected  surface  area,  see  [6]. 

We  suppose  that  the  plate  is  loaded  along  its  faces  3+R(e)  by  forces 
3 

e  (0,0, f (x;x/e))  per  unit  surface  area,  and  that  the  body  force  is 
2 

e  (0,0,F(x;jc/e))  per  unit  volume,  where 


(2.6)  f(x;n)  and  F(x;n)  are  bounded,  periodic  in 

n  ,  and  even  with  respect  to  • 

The  equations  of  elastostatic  equilibrium  for  the  clamped,  e-dependent, 
three-dimensional  plate  are 
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(2.7) 


(2.8) 


(2.9) 


-  3 


f 0  i=1’2 

”  '(A  1=3 

{0  1-1,2 

1-3 


u  =  0 


In  R(e) 


on  3+R(e) 


on  3qR(e)  . 


This  scaling  of  the  loads  ensures  that  ii  stays  bounded  as  e  -*  0 
Notice  that 


(2.10) 


c  e  e 

are  odd  ;  u^  is  even 


c  e  g 

0  , o „ _  are  odd  ;  o  ,  is  even 

ag  33  a3 


with  respect  to  ,  as  a  consequence  of  (2.6);  X£  will  denote  the  space 

of  all  admissible  displacements  that  obey  these  symmetries: 


(2.11) 


Xe  '  < 


—  3qR(e)  "  °  J  VU2  are 


u  €  H  (R(e))  : 


odd  and  is  even  in 


where  H  (R(e))  is  the  space  of  (vector-valued)  functions  with  square- 
integrable  first  derivatives. 

The  restriction  to  even  loads  is  merely  a  matter  of  technical  convenience. 
If  F  and  f  are  odd  In  »  then  the  solution  of  (2.7)-(2.9)  satisfies 


(2.12) 


II  e<ue)ll  2,e  *  c  e 


7/2 


In  case  |v^|  >,  c  >  0  on  3+R(e)  (i.e.  3+R(e)  has  no  vertical  parts), 
one  can  prove  (2.12)  by  taking  the  inner  product  of  uE  with  (2.7), 


integrating  by  parts,  and  using  a  Poincar£  inequality  on  each  vertical  line. 
The  proof  in  the  general  case  is  similar,  but  it  requires  the  methods  of 
section  4.  We  shall  show  here  that  for  even  loads 


(2.13) 


II  e(u£)  || 


3/2 


2,e 


wherever  the  "mean  load"  is  nonzero.  Since  the  problem  is  linear,  any  load 
can  be  decomposed  into  its  even  and  odd  components;  by  (2.12)  and  (2.13), 
the  even  part  is  the  one  that  produces  the  dominant  strain. 
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3.  Review  of  the  model 

The  model  presented  in  [6]  provides  the  initial  terms  of  an  asymptotic 
expansion  for  the  displacement  vector,  and  -  most  importantly  -  an  equation 
for  the  limiting  vertical  displacement  of  the  midplane.  This  equation  has 
in  each  case  the  form 


(3.1) 


9  o<M  O  *  3  Xw> 
aB  aByo  yo 


F  , 


where  F  =  F(F,f)(x)  is  the  rescaled  mean  vertical  load  (see  (5.1)  for  the 
precise  definition  of  F  when  a  =  1;  for  a  ^  1  see  [6]).  The  tensor 
M^Byd  represents  the  "effective  rigidity"  of  the  plate;  it  satisfies  the 
usual  symmetries 


m  =  M  =  m  =  M 
aBy<5  Bay6  aB<Sy  yfiaB  ' 


and  it  is  positive  definite  in  the  sense  that 


Ma8y6  ««3  *  c  I*1 


for  symmetric  tensors  E  „  .  The  formula  for  M  _  .  depends  on  whether 

’aB - aBy  6  — c - - - 

a  >  1  ,  a  =  1  ,  £r  a  <  1;  in  each  case  it  is  determined  by  h  through 
the  solution  of  certain  "cell  problems"  with  periodic  boundary  conditions. 


3A.  The  case  a  <  1. 


Let 


denote  the  set  of  functions  which  are 
periodic  with  period  L  =  (L^,L2)  ,  with  square  integrable  derivativesof 
order  $  2  .  The  auxiliary  functions  4>°^(n)  are  in  this  space,  and  are 
characterized  (modulo  a  constant)  by 


10 


and 


11 


(3.6)  lyd)  '  BljU  EU(i>  • 

Let  Q  denote  the  rescaled  period  cell  determined  by  h  , 


(3.7) 


Q  =  (n  ;  |nj  <  ~2  ,  |n3l  <  h(Q)}  ; 


and  let  F  denote  the  vector 


3/1  ' 

\  9  | 

<  l  > 

1  i_  .  . 

1-  K2 

B33y6 

82  /I  \\ 

3nx  \  211  a^B , 

r  n3  3n2  ’ 

(IVS/ 

I’TVb 

h  2^3 

B3333 

3n  3n-  V^vvi 

Y  e 

ctB  1 

The  auxiliary  functions  ^  €  H  (Q)  are  periodic  in  n  ,  and  they  satisfy 

(3.9)  =  -  | 

for  any  ^  €  H^Q)  which  is  periodic  in  n  .  The  tensor  is  given 

by 

■  #3iU«]  -  L"3SaSX»EXu'±YiX'5  ■ 


The  lowest  order  terms  in  the  displacement  vector  are 


u*  -  -  x,3  w  +  e2$aB(x/e)3  Qw 
Y  3  y  Y  —  oiB 


*  1,  .2  33a8  ,  .  2,aB  ,  . 

u3  '  "  +  I(,3>  Fj—  3a8”  +  *  *3  (i/=)3a6» 


(3.11) 
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3C.  The  case  a  >  1. 

We  define  a  tensor  »  ^or  use  this  section  only,  by 


Co363  "  C3o63  ’  Ce33B  "  C3a3B  "  0 


CljU  ‘  BijW  otherwise. 


For  any  function  £(n)  ,  £(&)  will  denote  the  associated  "stress", 


(3.12) 


hi  <*>  ■  C1,U  \£<±) 


ijk£ 


Let  Q  and  be  as  in  (3.7),  (3.8),  and  let  V  be  the  space  of 

2 

functions  j[^  €  L  (Q)  such  that 


is  periodic  in  n  ,  depends  only  on 


and 


|  Iij  (£)  Ei j  K  ®  • 


rtft 

When  a  >  1  ,  the  auxiliary  functions  are  in  V  ,  and  they  satisfy 


(3.13) 


Z1j(4ae>Eij(. -  -  j  ii;j(Ia8)E1;)(l)rfD 


for  every  €  y  .  The  tensor  giver»  by 


(3.14) 


M 


a&yS  L 


e1j<±tV5)*  . 


and  the  lowest  order  terms  in  the  displacement  vector  are 
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u*  =  -  x_3  w  +  e^+a$a®3  Qw 
Y  3  y  Y  aB 


(3.15) 


*  .  1,  ,2  B33aB  «  2,oB^ 

U3  =  w  +  2(x3)  —  3agw  +  c  3aew 


The  right  hand  side  of  (3.15)  is  evaluated  at  n  =  (x/ea  »x3/e)  . 

In  [6]  we  wrote  (3.13)-(3. lb)  in  a  sightly  different  form,  to  emphasize 

the  connection  with  homogenization  of  a  rough  surface.  The  functions 
08  08  33 

g  (n3)  ,  i|>  P(n)  ,  and  (n)  used  in  [6]  correspond  to  the  decomposition 


(3.16) 


(n3) 


aB 

g 


3 


2  B33ag 
B3333 


-  n3*“B(n) 


aB 


3n. 


(n3) 


33 

K  (n) 


One  can  characterize  <J;i^(*,n3)  as  the  solutions  of  certain  cell  problems 

ciB 

on  the  horizontal  slices  of  Q  ;  g  (^3)  may  be  expressed  in  terms  of 
certain  averages  of  . 

3D.  An  asymptotic  limit  of  the  case  a  ■  1. 

For  a  given  periodic  function  h(n>  ,  let  be  the  effective  rigidity 

of  the  associated  "a  <  1"  plate  defined  by  (3.3).  Let  denote  the 

aByo 

effective  rigidity  of  the  "a  *  1"  plate  with  thickness  variation 


h.(n)  *  h(n/A)  0  <  &  <  «  , 


i.e.,  is  as  defined  by  (3.10)  with  h  replaced  by  h.  .  We  show 

here  that  M*’*  as  4  -*■  ®  ,  if  h  is  smooth  enough;  the  proof  is 


similar  to  Nordgren's  convergence  argument  [9].  A  related  result  in  the 


14 


context  of  Laplace's  equation  may  be  found  in  [2] 
Proposition  3.1:  If  h  i  c  >  0  _is  smooth  enough 


(3.16) 


lim  M1’*  - 

^  agyfi  a6y6 


Proof; 

Let  denote  the  solution  of  (3.9)  with  1 

which  is  periodic  with  period  4L  .  We  introduce 

(3.17)  na  “  ^  na  ^a=1*2^  »  n3  = 


which  range  over  the  4-independent  cell 


(3.18) 


Q'  =  (V  ;  |n;|  <*|lq  ,  a  =  1,2  ; 


and  note  that 


EaB<*> 


24  ^  j  3n; 


(3.19) 


.  a*,  ,  a* 

t.  / ,  x  1  3,1  a 

Ea3(^  "  2l  3V  +  2  3^ 
a  j 


E33<^  “  TnT 


then 

i  replaced  by  h  , 
o 

the  rescaled  variables 
i3  » 

|n^|  <  h(g') ) 


Let  i|>aE(n')  be  the  solution  of  (3.2)  with  thickness  h(n')  ,  and  define 
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(3.20)  <£%’)  -  42(0,0,*aB)  +  4(-n’  ,  -  n’  ^  ,o) 


(n  n  h  '  \2  B33r<5  32/B 

V’ *2  3  *3333  3n;3n« 


We  shall  show  that 


(3.21) 


lim  [  =  0  , 

Jq'  4 


provided  that 


(3.22) 


*“%■> 


has  bounded  derivatives  of  order  5  4  , 


which  is  true  for  sufficiently  regular  h  .  Assertion  (3.16)  follows 
immediately  from  (3.21),  since 


(3.23) 


M 


1,4 

aBy6 


and 


lim 

A-*® 


L1L2 


Ma<1 

a6yfi 


Y 


For  fixed  a  and  3  ,  define  a  tensor  T^Cn1)  as  follows:  for 

6  €  (1,2}  , 


.2  ,a8 

t  -  5  _ ,  _  5  3  <p  , 

yS  oBy^I  y6po  3n’3n'  ^3  ’ 

P  a 


for  y  -  1,2, 
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\3  “  * 


1i 


7  7  *3,a6 

((^)2-h2(n')]B 

0  O  0 


r  i. 

aBy6  3n' 


2,aB 


-  B 


it 


v6pa  an'Sn' 
p  o 


3h 

3n 


;} 


and  is  the  solution  of 


(3.24) 


:_-33  .  -1  !!i3 

3n;  4  3n'  ’ 


331, 


no“+h(n') 


+  ."2  /_  =  .  3h_  3h_.  _  s  &VB  ,  3h  3h  " 

—  y  aBf<5  3tt^  3n^  Y<Spo  3rip3n^  3n^  3n^ 


One  verifies  the  consistency  condition  for  (3.24), 


£ 


h  *1x1  ,,  m  9-l  /=  .  3h_  3h_  ,  aVB  .  3h  3h 

-  j_h  3nj  °3  \  3r^  3n6  Y^PO  3np3no  8nY 


by  means  of  (3.2).  A  straightforward  computation  shows  that 


(3.25) 


lijCjfcf+JT)  -  ttJ  +ou_1)  , 


and  that 


<3-26)  aST  ’u  '  4'1  +  T5J  ■  0  ln  O' 

J  a  J 


l+JVh| 


11  1  \l+»  |vt 


'  0  °n  9iQ’  ' 


where  3+Q'  are  the  upper  and  lower  faces  of  Qf  ,  and  v_’  is  the  outward 
unit  vector  normal  to  3  Q*  . 

r1  a 
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4.  Integral  estimates 

This  section  establishes  certain  integral  inequalities  for 
u  €  H*(R(e))  .  We  consider  only  the  case  a  *»  1  ,  i.e. 

R(e)  *  lx  :  x  {  !1  ,  jx^j  <  eh(x/e)} 

(see,  however.  Remark  4.2  at  the  end  of  the  section).  Our  method  is  to 
decompose  R(e)  into  £?(e  )  subdomains,  each  with  diameter  of  order  c  , 

and  to  apply  Korn's  inequality  on  each  subdomain. 

We  begin  by  reviewing  Korn's  inequality  on  the  unit-sized  domain 

L 

Q  =  {x  :  |xa|  <  a/2  ,  |x3(  <  h(x) }  . 

R  is  the  space  of  rigid  motions, 

X  R  <=>  Y^OO  *  cijxj  +  ^i  *  ^or  some 

3 

ci  €  R  and  some  skew- symmetric 

matrix  c  . 


Vu  denotes  the  (nonsymmetric)  tensor  9^u^  »  and  e(u)  denotes  the 
(symmetric)  strain  tensor  y(3jUj+9^Uj)  * 


Lemma  4.1;  For  any  u  €  H  (Q)  there  exists  jr  €  R  such  that 


(4.1) 


[  iv(u-^)|2dx  $  C  [  |e(u)|2dx 

Jn  in 
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Proof :  This  follows,  for  any  Lipschitz  domain,  from  the  results  in  [5].  | 

Recall  that  R(e)  is  defined  by  (2.5).  For  each  pair  of  integers 
(k,-t)  let  R^£  *  R^(e)  denote  the  period  cell  centered  around 
(keL^.-deLj)  , 


eL, 


=  <x  :  (xj-keLjJ  <  ~~  ,  jx2~£cL2|  < 


eL 


_2 
2  ’ 


|x3|  <  eh(x/e)} 

Rescaling  (4.1)  and  (4.2)  yields  the  following  result. 

Lemma  4.2:  For  any  u  €  H  £qc(R(c))  and  any  pair  (k,£)  there  exists  y**'  € 
such  that 


(4.3) 


|v(u-x^)  |2dx  $  C  f  |  e  C u.)  | 2 cbc 

V  V 


and 


(4.4) 


f  (u-x^l2^  ?  Ce2  f  |e(u)[2cfx 

V  Jru 


The  constant  C  i£  (4.3),  (4.4)  depends  only  on  h  . 


Let 


k£,  .  k£  .  kl 
X  (x)  “  £  x  +  d  , 

\JL  3 

where  <1  €  F  and  £  is  a  skew-symmetric  matrix.  Clearly 


20 


(4 


c.  4,  k+l,f  k£,2  ^  2 1  ,  k+lf£  VZ  ( 2 

,5)  c  |c  -c  I  +  e  |?L  -d  I 


S  c  |  |^k+1*^(x)  -  2^(x)  |2dx2dx3 

S  C  j  (|lk+1'£-u|2  +  lu-jr^l2)^^  , 


where  the  integrals  are  over  the  interface  between  and  ^ 


eL„ 


{x  :  Xj  -  (k+%) eL^  ,  |x2-£eL2|  <  -y-  ,  jx3(  <  eh(x/e)} . 


One  has  the  trace  estimate  (on  any  Lipschitz  domain) 

(4.6)  f  | w| 2ct&  $  C  [  (|Vw|2+|w|2)dx 

J3Q  •'Q 

for  all  w  €  H^(Q)  .  Rescaling  (4.6),  and  combining  the  result  with  (4.5), 
we  obtain 

(4.7)  e4|ck+1,£-ck*|2  +  e2|dk+1,£-dk*|2 

$  C(ef  1  V(ju-xk+1,>t) \2dx  +  e-1  f  |u-xk+1,*|2clx 

^+1,*  ^c+l.i 

+  e  f  ^(u-x1^)  |2dx  +  e_1  [  lu-x^l2^)  . 

\e 


A  combination  of  (4.3),  (4.4)  and  (4.7)  gives 


2.  k+M  ke.2  .  ,,k+l,£  kd.i 2 

:  |£  -c_  I  +  |d_  ’  -d.  |  * 


Ce 


-1  l 


^+1,^  \z 


I e (“)  I  dx 


(4.8) 
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Similarly,  we  have 


,,  n.  2,  k,£+l  k£,2  -  .k,£+l  M,2 

(4.9)  e  |c  ’  -£  |  +|£  -d_  I 


<  Ce 


i 


Vf+l"  \t 


j e (u) (  dx 


Proposition  4.1:  For  any  £  €  H^RCe))  with  £  =  0  on  9»R(e)  , 


(4.10)  ||u||  2>£  +  ||Vu||  2j£  *  Ce"1  || e(u)  ||  2> 


Proof : 

Extend  u  to  R(e)  by  letting  it  be  zero  outside  R(e)  ,  and  let 

be  the  rigid  motions  introduced  in  Lemma  4.2;  notice  that  =  0 

if  R^,  fl  R(e)  =  0  .  Let  £(x)  and  6(x)  denote  the  piecewise  bilinear 

,  kl  ,  M 

interpolants  to  o_  and  d^  ,  i.e. 


CTij^keLl’^eL2^  =  Cij  ’ 


k l 


6i(keL1,£eL2)  *  d^  , 


and 


(x)  ,  <£(x)  are  bilinear  functions 


on 


(x  :|x1  -  (k4i)cL1|  <  ~  ,  |x2  -  (£+  |)eL2|  < 


for  each  pair  of  integers  (k ,t)  . 


It  is  standard  that 
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n  uj-1  p  \c(  ,2  i  ,k,£+l  ,k£. 1 2 . 

,  |«|2<ix  s  c  l  (li  '  -i  cl  +  \±'  -i  I  1 

Jr2  ~  *,l 

and  an  analogous  ine<u.lit,  holds  for  /|Vc|2&  .  It  follows,  using 


(4.8)~(4.9),  that 


M2d?  +  j^2  lyAi2^  *  Ce'1  ^  |e(u)j2dx  . 


Since  a  and  _6  are 


compactly  supported,  we  conclude  by  Poincard’s 


inequality  that 


j^2  lfi|2<k  +  r2  Ul2rf5  -  Ce  1  fR(e)  le(-)'2<^  * 


and  hence 


(4.11) 


T  [e2|cU|2  +  jd1^!2)  $  Ce"3  |e(u)|2ck 

jRfr) 


Since  Vj£ 


,k£  =  ck*  ,  (4.3)  may  be  rewritten 


[  lvu-cU|2(fa  5  C  f  |e(u)|2*t  , 

v  ■■  -  V 


which  leads  immediately  to 


(4.12) 


f  i V u  1 2 $  c  (f  |e(u)j2dx  +  l  e3|cU|  )  . 

)  R  (  c- 1  VR(e)  k>^ 


Si.il.rly,  since  |  S  C(|cM|+|d“|)  for  every  X  €  R(e)  ,  (*•«) 


leads  to 


(4.13)  |uj2dx  <  C  c2  [  |e(u)|2dx 

JR(e)  L  Jr(c) 


+  l  (|c“|2+|-k*|2)l 

k  ,1  J 


A  combination  of  (4.11),  (4.12)  and  (4.13)  gives 


R(e) 


(|u|2+|Vu|2)dx  <  Cc"2 


R(e) 


I  e (u) I  dx  , 


which  is  equivalent  to  (4.10). 

Recall  that  the  space  X^  is  defined  by  (2.11). 


Proposition  4.2:  For  any  €  X  , 


(4.14) 


l  "Uc+X33aU31!  2,t  *  Ce|le(^"  2 , e 

a=l 


k£ 

Proof:  When  u  €  X^  ,  one  may  choose  the  rigid  motions  of  Lemma  4 

to  have  the  same  symmetry  properties,  i.e. 


(4.15) 


,k-£  ,Vl  kl 
dl  ‘  d2  ‘  c12  ‘  ° 


for  each  k,£  .  By  (4.4)  and  (4.15), 


,  U  .2,  f  |  k£,  2  , 

u  -c  -x.  dx  =  u  -y  ax 

1  a  a3  3'  -  J  1  a  'a  1  - 


*U 


hi 


$  C  e2  [  |  e(ii)  |  2dx 

V 


By  (4.3),  on  the  other  hand. 
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Remark  4 . 1 

Inequalities  (4.10)  and  (4.14)  are  sharp  in  their  dependence  on  e  . 


For  (4.10),  one  sees  thi^  by  considering  u  =  (-x^a^w.-x^a^w.w) 
w  =  w(x1,x2)  .  For  (4.14),  one  uses  u  =  (0,0, w)  . 


where 


Remark  4.2 

The  estimates  in  this  section  may  be  generalized  considerably.  We 
assumed  that  u  vanishes  on  3gR(e)  to  simplify  its  extension  to  R(e)  . 
One  verifies,  with  a  little  more  effort,  that  Propositions  4.1  and  4.2 
remain  valid  without  this  condition  (modulo  a  rigid  motion,  in  the  case  of 
(4.10)).  The  argument  presented  here  also  works  in  the  case  a  <  1  ;  it 
applies,  moreover,  even  if  R(e)  is  only  a  Lipschitz  domain;  and  the 
periodicity  of  the  domain  is  not  essential. 

The  case  a  >  1  is  more  subtle;  we  do  not  know  nontrivial  conditions 
on  h  which  assure  (4.10)  or  (4.14)  for  that  scaling.  The  methods  of  [l] 
and  [4]  may  be  relevant  in  that  case.  | _ 


v 
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5 .  An  averaging  lemma 

Our  attention  remains  restricted  to  the  case  a  *  1  .  Q  denotes  the 
rescaled  period  cell  (3.7);  0  is  the  mean  thickness 

0  =  2M[h]  ; 

and  3+Q  is  the  "non-periodic"  part  of  3Q  , 

L 

3+Q  =  3Q  PI  (n  :  |n  j  <  -y-  ,  a  =  1,2}  . 

For  any  pair  of  functions  G(x;rj)  and  g(x;n)  which  are  periodic  in  n  , 
we  define 

(5.1)  F(G,g)  (x)  =  r—  1  I  Gdn  +  f  gdiin 

12  Uq  J3+Q 

(di  denotes  surface  measure).  Our  goal  is  the  following  result. 


Proposition  5.1:  Suppose  that  G  and  g  have  derivatives  in  x  of  order 
S  2  which  are  C^,a  and  C^,a  in  q  uniformly  in  x  ^  respectively.  Then 
for  any  u  f 


(5.2) 


|  G(x;x/e)u,dx  +  e  g(x;x/e)u_c& 

■'R(e)  *  J3+R(e) 


-  e 


-1  f  F(G, 

JR(e) 


g) (x)u.dx 


S  Ce3/2  ||e(u)  | 


2,c 


The  constant  C  depends  on  G  ,  g  ,  and  h  ,  but  not  on  e  . 

The  essence  of  Proposition  5.1  is  the  following:  if  F(G,g)  =  0  , 
€ 

and  if  w  solves 
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then 


(5.3) 


(°ij (wC)) 


r o  i-i,2 

lc2G  i-3 

CO  1-1,2 
\e3g  i-3 


e 


w 


0 


in  R(c) 


on  3+R ( e ) 
on  3QR(e)  , 


||e(we)  |j 


2,e 


S  Ce 


7/2 


Indeed,  if  G  and  g  are  even  in  ,  one  proves  (5.3)  by  substituting 

w  for  u  in  (5.2)  and  integrating  by  parts.  If  G  and  g  are  odd 

then  (5.3)  is  the  same  as  (2.12). 

Before  beginning  the  proof,  we  introduce  some  more  notation.  Given 
a  pair  G  ,  g  with  F(G,g)  =  0  ,  we  say  solves  the  cell  problem 
associated  to  G  and  g  "  if 

(5. A)  [  L  .  (!.)E< ,  0|0dn  -  I  G\Ji_c(n  +  [  gij>  dA(n) 

/q  J  /q  J  J  a+Q  J 

for  every  ^  €  K^(Q)  which  is  n  -  periodic. 

Recall  that  E($_)  and  £(£)  are  defined  by  (3.5)  and  (3.6).  One  verifies 

easily  that  (5.4)  has  an  r>  -  periodic  solution,  unique  up  to  a  translation. 

2  a 

Since  we  have  assumed  that  R(e)  is  a  C  ’  domain, 


(5.5) 


»!«<*>»  i,.  ‘  “II  c||  o.a+  Hsli  l.„> 


All  norms  are  on  the  rescaled  period  cell  Q  or  the  boundary  3+Q  .  The  constant 
C  depends  only  on  Q  and  not  on  jt  ,  which  occurs  in  (5.4)  —  (5.5)  as  a 
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parameter.  If  we  define 


(5.6) 


t..  -  l  (£)(x;x/e)  , 
J  ij 


then  (5.5)  leads  to 


(5.7) 


Proof  of  Proposition  5.1 

It  suffices  to  consider  the  case  F(G,g)  =  0  ;  the  general  case  will 
follow  by  considering  G'  =  G  -  0  ^F(G,g)  ,  g'  ■=  g  . 

Let  ^  solve  the  cell  problem  (5.4)  associated  to  G  and  g  ;  let 
t  be  as  in  (5.6);  and  let  u  €  .  By  Green's  formula, 

(5.8)  f  Gu-dx  +  c  gu.d4 

JR(e)  J3+R(e) 


=  c 


'R(e) 


3x 


8  iB 


Ujdx  +  E 


R(e) 


Tij  eij(-d- 


where 


tion. 


_3 

3x 


to 


3G_ 

3x 


I  -ir  I 

B  IB  B, 
Notice  that  r — 

. 


and 


,0,a 


Y 

and 


3 


4, a 


*Y 


in 


(£)  is  evaluated  at  n  *  .x/e  after  differentia- 
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■  £  is  the  stress  of  the  cell  problem  associated 
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since  the  x  -  derivatives  of  G  and  g  are  assumed 
,  uniformly  in  x  , 


(5.9) 


l-sr  I  I 

3XB  iB 


2,e 


5  Ce 


1/2 


We  estimate  the  various  terms  in  (5.8)  separately.  First, 


(5.10) 


5  Ce1^2  ||e(u) 


2,e 


by  (5.7)  and  Haider's  Inequality.  Next, 


(5.11) 


- —  J  Q*(u  +x  3  u. ) dx  $  Ce  u  +x_3  u 
n/  .  3x~  *-aB  a  3  a  3  —I  11  a  3  a  3"  2,e 

R(e)  8  ' 


*  Cc3/2  ||  e(u)  II  2 1  E 


by  (5.9),  Holder's  inequality,  and  Proposition  4.2.  Finally, 


(5.12) 


[  IT —  I  bxt8 

JR(c)  3xB  “B  3  «  3  - 


s  Ce3/2  ||Vu|| 

*-  y 


<  Ce1/2  || e (u)  ||  2^e 


by  (5.9),  Holder's  inequality,  and  Proposition  4.1.  Combining  (5 . 10) -(5 . 12 ) , 
we  conclude  that  the  right  side  of  (5.8)  equals 


(5.13) 


f 


E  I  £  u3d-  +  <?(e3/2  IMh)II  ?  c> 

JR(e)  3XB  36 


2  ,e 


At  this  point  we  need  the  following  identity,  which  will  be  proved  later 


(5.14) 


dri  =  0  . 


Q  B  38 


This  means  that  f(— £  ,  0  )  =  0  .  Repeating  the  above  argument  with 

\  XB  36  /  „ 


B  36 

the  cell  problem  associated  to  ^ 


l  and  0  ,  we  conclude  that 
6  3B 


(5.15) 


|  JL  l  u3dx  -  c2  | 
JR(e)  dx6  3B  J  JR(e) 


H  u-dx 
3Xp  3B  3  - 


+  0(c5/2||e(u)||2>e)  , 


where 


are  the  stresses  associated  to  the  new  cell  problem.  (We  use 
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here  our  hypothesis  on  the  second  derivatives  of  G  and  g.)  By  Holder's 
inequality.  Proposition  4.1,  and  the  analogue  of  (5.9), 

(5.1.6)  -3  u  dx  $  Ce  1/2  ||e(u)|^  . 

R(e)  3  0  36  3  2,c 

Combining  (5.8),  (5.13),  (5.15),  and  (5.16),  we  obtain  (5.2). 

It  remains  to  prove  (5.14).  Substituting  ^  =  (n^.O.O)  into  (5.4) 

gives 

|  I3l^)dD  ■  |  ^ij(i)Eij(^)dil 
=  0  . 

Since  Q  is  independent  of  x  ,  it  follows  that 

fQ  ^31dD  “  8^  fQ  ^31dQ  =  °  * 

9  v 

The  corresponding  assertion  for  follows  using  £  *  (0,n3,0)  , 

and  summation  leads  to  (5.14).  I 
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6.  Convergence 
Let  w  solve 


VWV>  "  F("'f>  ln  ° 

w  =  3  w  =  0  on  3!1  , 

n  ’ 

where  M  .  .  Is  defined  by  (3.10);  let  u*  be  as  in  (3.11);  and  let 
o8y6  — 

ue  solve  the  three  dimensional  elasticity  problem  (2.7)-(2.9).  We  shall 
prove  that  u*-uE  converges  to  zero  in  energy,  and  that  w  is  really  the 
limiting  vertical  displacement.  In  addition  to  the  regularity  hypotheses 
on  0  and  R(e)  ,  formulated  in  section  2,  we  must  assume  that 

(6.1)  All  x-derivatives  of  order  £  2  of  F  and  f 
are  C°’a  and  C'*''0  in  n  ,  respectively 
(uniformly  in  x  ). 

By  [6],  M  is  positive  definite;  it  follows  that 

apyo 

(6.2)  w  has  bounded  x-derivatives  of  order  $  6  . 

Let  f;(t)  ^  C^(0,®)  with  ?(t)  =  0  for  t  $  1/2  and  s(t)  =  1  for 
t^l;  we  define 

£^(x)  =  c(e  1dist(x,3ft))  for  x  €  Q  , 

and 

(6.3)  t/  -  (-x^w.-x^w.w) 

g 

+  <O.O.j(x3)2  . 
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where  Is  as  in  (3.9).  Notice  that  u^  €  X^_ 

2 

Since  R (c)  is  assumed  to  be  a  C  '  domain,  a  standard  regularity 
5  2 

result  shows  that  <j>_'  (  C  1  (Q)  .  In  particular  the  functions 


V±YS>  -  > 


are  and  C^’a  ,  respectively.  We  shall  use  this  fact  repeatedly, 

sometimes  without  direct  mention,  in  what  follows. 


Lemma  6.1:  The  functions  u^  ,  u*  satisfy 


(6.4) 


|e(u*-u*)||  £  C c2 


Proof: 


Since  and  ^  (^  p)  are  bounded  functions. 


I  v^c  I  Ci£/e>  li  2,  e  S  c 


||  (l-Ce)iaP(x/e)  ||2  £  S  C e 


||  (l'4G)E(iae)(x/e)||2  e  s  Ce 


The  estimate  (6.4)  follows  easily,  using  (6.2). 


Lemma  6.2:  For  each  &  ,  y  ,  6 


(6.5) 


|  [3g(iY6)rfD  “  0 


Proof : 


One  argues  as  in  the  proof  of  (5.14),  using  (3.9)  instead  of  (5.4)  | _ | 
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on  the  other  hand 


y&\  . 


(6.9) 


i  _  v($.  afl33  r  (*,7°) 

B  E  Y  )  *  laR^~  ^  ~  ^33  ^ 

afUu  Xp  ~  ap  3333 


combination  of  (6.8)  and  (6.9)  yields  (6.7). 

One  easily  verifies  that  (3.9)  is  eqnivalent  to 


JL  y  (<bY6)  =  0  in  Q 
"j  - 


■  n3Bi6y6v6  0n  *+' 


on  3  ,Q 


jf  t  is  defined  by 


...  - 


(6.10)  Ttj  "  1  X3  ijy«  ^ 


one  computes  that 

(6.11)  3jTij  =  l_X3"iSy6 

in  R(e)  ,  and 

(6.12) 


on  3_jR(e)  ■ 


a.t..  -  I16(t,!><i'c,,38rd'‘ 


■  ° 


Proposition  6 For  anx  v  e  xc  ’ 

,  f  T  e  (v)dx  =  e2©"1  [  F<F’f)  v3d* 

(6-13)  ]  }  Tij  ij  -  -  jR(e) 
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+  0(e5/2  lie(v)||  2>e>  • 
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Proof : 

By  Green's  formula  and  (6.12), 


T..e..(v)dx  =  -  3 . (t.  ,  )v  dx 

J»(t)  ij  -  ~  W>  5  13  1  - 


We  apply  Proposition  5.1,  using  (6.11),  (6.5),  and  the  fact  that  B_  =  0, 

Jtsyo 


to  see  that 


}  MT3i)v3^  *  Cc5/2  He(^ll  2  , 

' R(e)  3  3  ’ 


Writing  v  =  -  x,3  v,  +  (v  +x-3  v-)  ,  and  applying  Proposition  4.2,  we 
a  JaJ  a  j  a  J 

obtain 


R  (e) 


3  .  (r  , )v  dx  =  - 
J  aj  a  — 


R(c) 


x_3  (t  .  )3  v.dx 
3  j  aj  a  3  ~ 


+  0(eV2  ||e(v)||  2>e)  . 


A  combination  of  these  results  yields 


(6.14) 


L(e) 


Tijeij(v)<4  = 


R(e) 


x  3  (x  , )  3  v  dx 
3  j  aj  a  3  - 


+  0(t5/2  ||e(v)||  2>e) 


We  use  Green's  formula  again: 


(6.15) 


x,3. (t  . )  3  v~dx  =  - 
3  j  aj  a  3  - 


x  3  (t  ,)v  dx 
3  aj  aj  3  — 


R(e) 


x_3  (t  )vev_d6 
3  j  aj  a  3 
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Now, 


x.,3.(t  .)vE  =  e2(l+ll)  ,  with 
3  j  aj  a 


1  n3BaeY6Va|n=x/E36Y6W  ’ 


11  =  n3  U(lY6)V^ 


l=x/ E 


3Sy6W  ’ 


and 


-  x.9  ,(t  )  =  c III  +  e  IV  ,  with 

3  aj  aj 


111  =  -  *3  W  W?  > 

a 


T^X/  E 


aeY«"  • 


rv  '  "'jW'Y  Lb(1y6,) 


n=x/  e 


3aBy6W 


By  Lemma  6.3  and  Proposition  5.1, 


(6.16)  e 


R(e) 


(Ill)v^dx  +  z 


(I+lDv^di  =  0( e^2  1 1  e ( v)  1 1  ^  E ) 


3+R<e) 


similarly,  by  Lemma  6.4,  Proposition  5.1  and  (3.10), 


(6.17)  e 


R(e) 


(IV)v^dx  =  e  6 


2  -1 


Ma6Y63a8y^  ^dx  +  0(e7/2  ||  e  (v)  ||  ^ 


R(e) 


Since  M  „  r  is  constant  and  w  satisfies  (3.1),  (6. 15)  —  (6. 17)  imply 
oipyo 


(6.18) 


R(c) 


3,(i  ,)8  v.dx=  e20-1  f  F(F,f)v.dx  +  0(c5/2  ||e(v)||  _  )  . 

j  aj  a  j  Jr(£)  5  z*c 


A  combination  of  (6.14)  and  (6.18)  yields  (6.13).  | 

We  are  ready  to  prove  the  main  result  of  this  paper. 

Theorem  6.1:  The  ansatz  tj  ,  defined  by  (3.11) ,  and  the  displacement  uc 
defined  by  (2.7)-(2.9),  satisfy 

(6.19)  ||e(u*-u£)  ||  <  Cc2  . 

^ 

Proof : 

We  shall  prove 

(6.20)  lle()i//-liC)||  7  i  Cz2  , 

z ,  e 

with  /  as  in  (6.3).  The  estimate  (6.19)  is  an  immediate  consequence, 
using  Lemma  6.1. 

To  prove  (6.20),  we  first  observe  that 

(6.21)  I|o(/)-t|U  S  Cc2 

^ ,  e 

where  x  is  defined  by  (6.10).  Indeed,  a  simple  computation  gives  that 

||o(u*)-T||  <  ce5/2  , 

^  ,  c. 

while  by  Lemma  6.1, 

|jo(u*)-a(u^)  ||  <  Cc2  ; 

(6.21)  follows  by  means  of  the  triangle  inequality. 


By  Proposition  6.1  and  (6.21) 
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(6.22)  o  .  . (u  )e . . (v)dx  =  c“ 

-  ij  ~  - 


-1  ' 

■  R(c) 


F-v3dx  +  0(c  || e (v)  ||  2  £) 


for  any  y  €  .  Also, 


(6.23) 


c  2  3 

o. . (u  )e. . (y)dx  =  c  Fv.dx  +  £  fv^di 

He)  1J  13  R(t)  J3+R(e) 


=  e20  1  [  Fv  dx  +  OCc7^2  ||e(v_>(|  ) 


by  (2.7)-(2.9),  Green's  formula  and  Proposition  5.1.  Taking  v  =  u' -uE  , 
and  subtracting  (6.23)  from  fb.22),  we  conclude  that 


(u  -uC)e^.(u  -u  )dx  $  Ce  j|e(u  -uE)(|  2  ^  , 


from  which  (6.20)  follows. 


Remark  6.1 


Had  we  specified  the  e-dependent  boundary  condition 


3RQ(e)  |3R0(c) 


instead  of  (2.9)  then  the  introduction  of  u  would  not  have  been  necessary. 
The  above  argument  yields  j|e(u  -uE)||  $  Ce"*^2  when  u6  is  defined 

this  way.  |  | 

^  "i  /  *) 

One  verifies  readily  that  ||e(u  )  ||  9  ~e  “  whenever  F  ^  0  .  It 
follows,  using  (6.19),  that 


(6.24) 


C_1e3/2  f  ||e(uE)||  $  Ce3/2  , 

Z.  9\, 


with  C  depending  on  F  but  not  on  t  .  A  combination  of  (6.19)  and 
(6.24)  yields  the  relative  error  estimate 
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l!e(u*-uC] 
||e(uC)  || 


2 J  ..  1/2 
— <  Cc 


2,c 


A  similar  argument,  using  (4.10),  shows  that 


laVlL, 


c.« 

Iucl!2.c  ll^ll2>c 


One  may  also  compare  u^  and  w  directly: 


Corollary  6.1:  If  one  defines 


wC(x)  =  1 


2eh(x/ e) 


(-+eli(x/t ) 
-c.h(x/e) 


u3(x)dx3 


wherever  h(x/e)  /  0  ,  then 


(6.25) 


(f  |w-we|2h(x/e)cfx)1/2  ?  Ce1/2 
Jo 


Proof : 

We  consider  w 


e 

w 


to  be  defined  on  R(e) 


By  Poincar£' 


re  h(x/  e) 
'-eh(x/e) 


|wE-u3|2dx3 


s  Ce 


eh(x/e) 

-eh(x/e) 


3 

3x. 


uE|2dx3 


re h(x/e) 

■*  -eh(x/e) 


le33(HC)  |2^3 


inequality 


integration  over  ii  yields 
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I  c  e  1 2  .  2 

[uy-w  [  ax  f  Cc 


[ 

JR(e) 


c)  1 2cfx  *  Cc5  . 


One  computes  directly  from  (6.3)  that 


[  lw-ut  |  £  Cc3  . 

JR(e) 

Combining  these  two  estimates  with  (4.10)  and  (6.20),  we  conclude  that 
(f  |u-wC|2dx)1/2  <  Ce5/2  +  ||u^-u||!  ,  .  *  Cc  . 

Jr(c)  3  3  "’C 


It  follows  that 


|w-wE  |  2h(x/c)dx  =  ~~  [  |w-wE|Zc(x  ?  Cc 

JR(e) 


Remark  6.2:  If  h(n)  ^  c  >  0  ,  i.e.  if  the  plate  has  no  holes,  then  (6.25) 


becomes 


(  |w-wG| 
Jq 


2dx)1/2  *  cc1/2 
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